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Abstract 



The ultimate goal for future neutrino facilities is the determination of CP violation 
in neutrino oscillations. Besides |{7(e3)| 7^ 0, this will require precision experiments 
with a very intense neutrino source and energy control. With this objective in mind, 
the creation of monochromatic neutrino beams from the electron capture decay of 
boosted ions by the SPS of CERN has been proposed. We discuss the capabilities 
of such a facility as a function of the energy of the boost and the baseline for the 
detector. We compare the physics potential for two different configurations: I) 
7 = 90 and 7 = 195 (maximum achievable at present SPS) to Frejus; II) 7 = 195 
and 7 = 440 (maximum achievable at upgraded SPS) to Canfranc. We conclude 
that the SPS upgrade to 1000 GeV is important to reach a better sensitivity to CP 
violation iff it is accompanied by a longer baseline. In both Setups, the gain in the 
CP violation sensitivity with a previous knowledge of |{7(e3)| is apparent. 



1 Introduction 



Neutrinos do have masses and mixings. Present evidence [1,2] from neutrino osci- 
llations is consistently interpreted in terms of two independent mass-differences and 
mixings: the so-called atmospheric sector (2,3) and the solar sector (1,2). The initial 
discovery of the zenith effect with atmospheric neutrinos led to the (2,3) sector in 
neutrino oscillations, later confirmed by long baseline experiments with accelerator 
neutrinos. The solar neutrino solution to the historical solar neutrino problem led 
to the (1,2) sector from neutrino oscillations in solar matter, later confirmed by 
long-baseline vacuum oscillations with reactor neutrinos. Whereas the sign of Ara^ 
is thus known, the determination of the sign of Am^ needs the incorporation of 
matter effects and then |C(e3)| 7^ 0. The two mixings angles are large: 623 could 
even be 45°, whereas #12 is large although not maximal. 

The third connecting mixing |C/(e3)| is bounded as #13 < 10° from the CHOOZ 
reactor experiment [3]. The angle #13 remains thus undetermined. The approved 
reactor experiments Double CHOOZ [4] and Daya-Bay [5], as well as the second gen- 
eration of long-baseline superbeam experiments T2K [6] and NOVA [7] will address 
this point. A number of experimental facilities to significantly improve on present 
sensitivity and look for CP-violation have been discussed in the literature: neutrino 
factories (neutrino beams from boosted-muon decays) [8-10], superbeams (very in- 
tense conventional neutrino beams) [11—14], improved reactor experiments [15] and 
/3-bcams [16]. The original standard scenario for beta beams with low 7 = 60/100 
and short baseline L = 130 Km from CERN to Frejus with 6 He and l8 Ne ions 
can have a variant by using an electron capture facility for monochromatic neutrino 
beams [17]. New proposals also include the high Q value s Li and s Be isotopes in 
a 7 = 100 facility [18]. For the standard beta beam facility, a study of the physics 
reach as function of the boost and the baseline has been made [19]. In this paper 
we discuss the physics reach that a high energy facility for EC beams may provide 
with the expected SPS upgrade at CERN. In Section 2 we discuss the virtues of the 
suppressed oscillation channel {y e — > v^) in order to have access to the parameters 
013 and 5. To disentangle the CP-violating phase we emphasize the method of using 
energy dependence, as obtainable in the EC facility. In Section 3 we present new 
results on the comparison between (low energies, short baseline) and (high energies, 
long baseline) configurations for an EC facility with a single ion. Section 4 gives our 
conclusions. 



2 CP-even and CP-odd terms 

The observation of CP violation needs an experiment in which the emergence of 
another neutrino flavour is detected rather than the deficiency of the original flavour 
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of the neutrinos. At the same time, the interference needed to generate CP-violating 
observables can be enhanced if both the atmospheric and solar components have a 
similar magnitude. This can happen in the suppressed v e — > v„ transition. The 
appearance probability P{y e — > Vu) as a function of the distance between source and 
detector (L) is given by [20] 



where J = C13 sin 2#i2 sin 2#23 sin 2^13. The three terms of Eq. ([I]) correspond, res- 
pectively, to contributions from the atmospheric and solar sectors and their interfer- 
ence. As seen, the CP violating contribution has to include all mixings and neutrino 
mass differences to become observable. The four measured parameters (Am^ 2 ,^i2) 
and (Am23,#23) have been fixed throughout this paper to their mean values [21]. 

Neutrino oscillation phenomena are energy dependent (see Fig(T]) for a fixed 
distance between source and detector, and the observation of this energy dependence 
would disentangle the two important parameters: whereas |£/(e3)| gives the strength 
of the appearance probability, the CP phase 5 acts as a phase-shift in the interference 
pattern. In fact, a general theorem [22] states that, under the assumptions of CPT 
invariance and absence of absorptive parts, the CP-odd probability is odd in time 
or, equivalently, odd in L. As vacuum oscillations depend on L/E, this result 
implies that, for fixed L, the CP-odd probability is odd in E, whereas the CP-even 
terms are even in E. This is satisfied by Eq. ([Tj), with three contributions, the 
atmospheric, the solar and their CP-even interference, which are even functions of 
E, and the CP-odd interference, which is odd in E. These properties suggest the 
consideration of a facility able to study the detailed energy dependence by means 
of fine tuning of a boosted monochromatic neutrino beam. In an electron capture 
facility the neutrino energy is dictated by the chosen boost of the ion source and 
the neutrino beam luminosity is concentrated at a single known energy which may 
be chosen at will for the values in which the sensitivity for the (#13, 5) parameters is 
higher. This is in contrast to beams with a continuous spectrum, where the intensity 
is shared between sensitive and non sensitive regions. Furthermore, the known 
definite energy would help in the control of both the systematics and the detector 
background. In the beams with a continuous spectrum, the neutrino energy has to 
be reconstructed in the detector. In water-Cerenkov detectors, this reconstruction is 
made from supposed quasielastic events by measuring both the energy and direction 
of the charged lepton. This procedure suffers from non-quasielastic background, 





2 



from kinematic deviations due to the nuclear Fermi momentum and from dynamical 
suppression due to exclusion effects [23]. 

6 13 = 5 cleg 
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Figure 1: The appearance probability P(y e — > v^) for neutrino oscillations as a 
function of the LAB energy E/L, with fixed connecting mixing. The three curves 
refer to different values of the CP violating phase 5. The vertical lines are the 
energies of our simulation study in the EC facility. 

The above discussion proves that the study of neutrino oscillations in terms 
of neutrino energy will be able to separate out the CP phase 5 from the mixing 
parameters. A control of this energy may be obtained from the choice of the boost 
in the EC facility with a single ion. In order for this concept to become operational, 
it is necessary to combine it with the recent discovery of nuclei far from the stability 
line, having super allowed spin-isospin transitions to a giant Gamow- Teller resonance 
kinematically accessible [24]. Thus the rare-earth nuclei above 146 Gd have a small 
enough half-life to allow electron capture processes in the decay ring. This is in 
contrast with the proposal of EC beams with fully stripped long-lived ions [25]. We 
discuss the option of short-lived ions [17]. 
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3 Physics reach at different energies and baselines 



Electron Capture is the process in which an atomic electron is captured by a proton 
of the nucleus leading to a nuclear state of the same mass number A, replacing 
the proton by a neutron, and a neutrino. Its probability amplitude is proportional 
to the atomic wavefunction at the origin, so that it becomes competitive with the 
nuclear (3 + decay at high atomic number Z . Kinematically, it is a two body decay 
of the atomic ion into a nucleus and the neutrino, so that the neutrino energy is well 
defined and given by the difference between the initial and final atomic masses minus 
the excitation energy of the final nuclear state. In general, the high Z nuclear beta- 
plus decay and electron-capture (EC) transitions are very "forbidden", i.e., 
disfavoured, because the energetic window open in these channels does not contain 
the important Gamow- Teller strength excitation seen in (p,n) reactions. There are 
a few cases, however, where the Gamow- Teller resonance can be populated having 
the occasion of a strong "allowed" transition. For the rare-earth nuclei above 14e Gd, 
the filling of the intruder level h\\/2 for protons opens the possibility of a spin- 
isospin transition to the allowed level /19/2 for neutrons, leading to a fast decay. Our 
studies for neutrino beam capabilities have used the 150 Dy ion with half life of 7.2 
min, a Branching Ratio to neutrino channels of 64% (fully by EC) and neutrino 
energy of 1.4 MeV in the CM. frame, as obtained from its decay to the single giant 
Gamow- Teller resonance in the daugther 150 T6*. 

The parent radioactive ion is accelerated and then accumulated and storaged. A 
neutrino of energy Eq in C. M. that emerges from the decay in these conditions will 
be boosted in energy. This LAB energy is a function of the angle (9) of neutrino 
detection and Lorentz gamma (7) of the ion at the moment of decay and it can 
be expressed as E = Eg/ [7(1 — j3cos6)\. The angle 6 here expresses the deviation 
between the actual neutrino detection and the ideal detector position in the prolon- 
gation of one of the long straight sections of the decay ring. The neutrinos emerging 
from a boosted ion beam decaying by EC are concentrated inside a narrow cone 
around the forward direction. If the ions are kept in the decay ring longer than the 
half-life, the energy distribution of the Neutrino Flux arriving to the detector in the 
forward direction, in absence of neutrino oscillations, is given by the Master Formula 



d 2 N u 1 d 2 T u 

dSdE TdSdE u 
T V N U 



r -kL 2 



j 2 5(E-2 7 E ), (2) 



with a dilation factor 7 >> 1. It is remarkable that the result is given only in terms 
of the branching ratio for electron capture and the neutrino energy and independent 
of nuclear models. In Eq. ([2]), Ni ons is the total number of ions decaying to neutrinos. 
At the first oscillation maximum, with E/L fixed, Eq. ([2]) says that lower neutrino 
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energies Eq in the proper frame give higher neutrino fluxes. The number of events 
will increase with higher neutrino energies as the cross section increases with energy. 
To conclude, in the forward direction the neutrino energy is fixed by the boost 
E = 2"fEo, with the entire neutrino flux concentrated at this energy. As a result, 
such a facility will measure the neutrino oscillation parameters as a function of 
energy by changing the 7's of the decay ring (energy dependent measurement) and 
there is no need of energy reconstruction in the detector. In this situation, the 
experiment becomes a counting-rate of events. 
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Figure 2: Setup I as explained in the text. Fit for (#13, d) from statistical distribution 
of events with assumed values of the parameters. 
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Figure 3: The same as Figure 2 for Setup II as explained in the text. 
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For the study of the physics reach associated with such a facility, we combine 
two different energies for the same 150 Dy ion in each of two Setups. In all cases 
we consider 10 18 decaying ions/year, a water Cerenkov Detector with fiducial mass 
of 440 Kton and both appearance (vn) and disappearance (u e ) events. Setup I is 
associated with a five year run at 7 = 90 (close to the minimum energy to avoid 
atmospheric neutrino background) plus a five year run at 7 = 195 (the maximum 
energy achievable at present SPS), with a baseline L = 130 Km from CERN to 
Frejus. The results for Setup I are going to be compared with those for Setup II, 
associated with a five year run at 7 = 195 plus a five year run at 7 = 440 (the 
maximum achievable at the upgraded SPS with proton energy of 1000 GeV), with 
a baseline L = 650 Km from CERN to Canfranc. As explained in Section 2, the 
virtues of having at least two energies in a given Setup are that the two oscillation 
parameters #13 and S can be separated out. 
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Figure 4: Setup II. Sensitivity to #13 7^ for a two parameter (#13, 5) fit to the 
statistical distribution of events. 

For the Setup I we generate the statistical distribution of events from assumed 
values of #13 and 5. The corresponding fit with two parameters is shown in Fig. [2] for 
selected values of $13 from 8° to 1° and covering a few values of the CP phase 6. As 
observed, the principle of an energy dependent measurement (illustrated here with 
two energies) is working to separate out the two parameters. With this configuration, 
the precision obtainable for the mixing is much better than that for the CP phase. 
As seen, even mixings of 1° are still distinguishable. We emphasize that these results 
are obtained with a two-parameter fit, i.e., assuming that both (6*13, 5) are unknown 
quantities. 

At the time of the operation of this proposed Facility, it could happen that the 
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connecting mixing #13 is already known from the approved experiments for second 
generation neutrino oscillations, like Double CHOOZ, Daya-Bay, T2K and NOVA. 
To illustrate the gain obtainable in the sensitivity to discover CP violation from the 
previous knowledge of #13, we have reanalyzed the statistical distribution of events 
with the assumption of #13 already known in advance. In general, the precision to 
obtain 5 is then much better than that of Fig. [2] and the corresponding sensitivity 
for a CP-violation discovery is discussed later. 
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Figure 5: CP violation exclusion plot at 99% CL, if #13 is still unknown, for the two 
reference Setups: I (broken blue line) and II (continuous red line). 
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Figure 6: CP violation exclusion plot at 99% CL, if #13 is previously known, for the 
two reference Setups: I (broken blue line) and II (continuous red line). 

In the case of Setup II the longer baseline for 7 = 195 leads to a value of E/L 
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well inside the second oscillation (see Fig. [I]). In that case the associated strip in 
the (#13, 5) plane has a more pronounced curvature, so that the two parameters 
can be better disantangled. The statistical distribution generated for some assumed 
values of (#13, 5) has been fitted and the x 2 values obtained. The results are given in 
Fig. [3] for a two-parameter fit. Qualitatively, one notices that the precision reachable 
for the CP phase is better than that in the case of Setup I. One should emphasize 
that this improvement in the CP phase determination has been obtained with the 
neutrino channel only, using two appropriate different energies. One may discuss in 
this Setup II the sensitivity to discover # 13 7^ by giving the \ 2 fit, f° r each # 13 , to 
the value #13 = 0. This is given in Fig. [H Although it is somewhat dependent on 
the 5-value, we see that values of #13 > 1° are in general distinguishable for zero. 

The corresponding exclusion plots for CP violation in the two Setups are com- 
pared when both #13 and 5 are unknown. The sensitivity to discover CP violation 
has been studied by obtaining the x 2 fit for 5 = 0, 180° if the assumed value is 6. 
For 99% CL, the sensitivities to see CP violation in both Setups are compared in 
Fig. EJ In both cases, we assume a two-parameter fit, i.e., #13 previously unknown. 
For Setup I, a non-vanishing CP violation becomes significant for #13 > 4°, with 
values of the phase 5 around 30° or larger to be distinguished from zero. For Setup 
II, the sensitivity to CP violation is better and significant even at 1° in certain cases, 
depending on the hemisphere for the value of the phase d. 
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Figure 7: Setup II. CPV sensitivity for the statistical distribution of events depen- 
ding on a single parameter 5, assuming previous information on the value of #13. 

If #13 is previously known, the corresponding analysis for the sensitivity to dis- 
cover CP violation is presented in Figure [6l In this case, the x 2 fit is made with 
the single parameter 5. One may notice that the improvement in this sensitivity is 
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impressive, suggesting that going step by step in the determination of the neutrino 
oscillation parameters by means of several generation experiments is very rewarding. 
As in Fig. [5j Setup II provides better sensitivity to the discovery of CP violation 
than Setup I. In the best case, i.e., #13 already known at the time of the proposed 
experiment with Setup II, we give in Fig. [7] the sensitivity to discover CP violation 
for different \ 2 to be distinguished from 5 = 0, 180°. The result is so good that it 
enters into the regime of a precision experiment. 



4 Conclusions and Outlook 

The simulations of the physics output for an EC neutrino beam at different energies 
indicate: 

1) The principle of energy dependence to separate out the CP-even and CP-odd 
contributions to the neutrino oscillation probability works. 

2) The upgrade to higher energy in the SPS boost (E p = 1000 GeV) helps to 
have a better sensitivity to CP violation, which is the main objective of the next 
generation neutrino oscillation experiments, iff accompanied by a longer baseline. 

3) The best E/L in order to have a higher sensitivity to the mixing |C/(e3)| is 
not the same than that for the CP phase 5. Like the phase-shifts in interference 
phenomena, the presence of 5 is easier to observe when the energy of the neutrino 
beams enters into the region of the second oscillation. The mixing is better seen 
around the first oscillation maximum, instead. 

4) The previous knowledge on the connecting mixing #13 would greatly improve 
the sensitivity to CP violation discovery in this facility. This statement is valid in 
both experimental Setups: I of lower energy, shorter baseline, or II of higher energy, 
longer baseline. 

5) In the best configuration, i.e., with #13 known in advance and Setup II, the 
CP-violation sensitivity is of a few degrees in 5 for #13 > 1°. 

Besides the feasibility studies for the machine, most important for physics is the 
study of the optimal configuration by combining low energy with high energy neu- 
trino beams, short baseline with long baseline and/or EC monochromatic neutrinos 
with G He j3~ antineutrinos. 

Among the possible systematics associated with the proposed experiments, one 
should define a program to determine independently the relevant cross sections of 
electron and muon neutrinos and antineutrinos with water in the relevant energy 
region from several hundreds of MeV's to 1 GeV or so. 

The result of the synergy of Neutrino Physics with Nuclear Physics (EURISOL) 
and LHC Physics (SPS upgrade) for the Facility at CERN could be completed with 
the synergy with Astroparticle Physics for the Detector, which could be common to 
neutrino oscillation studies with terrestrial beams, atmospheric neutrinos (sensitive 
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to the neutrino mass hierarchy through matter effects [26]), Supernova neutrinos 
and Proton decay. 

The analysis shown in this paper indicates that the proposals discussed here merit 
RSzD studies in the immediate future for all their ingredients: Facility, Detector and 
Physics. 
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